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Abstract 

We report experiments investigating the behaviour of micron-scale 
fluid droplets jetted onto surfaces patterned with lyophobic and lyophilic 
stripes. The final droplet shape is shown to depend on the droplet size 
relative to that of the stripes. In particular when the droplet radius is 
of the same order as the stripe width, the final shape is determined by 
the dynamic evolution of the drop and shows a sensitive dependence 
on the initial droplet position and velocity. Lattice Boltzmann numer- 
ical solutions of the dynamical equations of motion of the drop provide 
a close quantitative match to the experimental results. This proves 
helpful in interpreting the data and allows for accurate prediction of 
fluid droplet behaviour for a wide range of surfaces. 

1 Introduction 

This paper presents a combined experimental and numerical investigation of 
the behaviour of micron-scale fluid droplets jetted onto chemically patterned 
surfaces!^ It was motivated by questions of relevance to ink-jet printing 
where substrates with chemical or physical defects can cause the expected 
spherical shape of jetted droplets to become distorted, thus affecting the 
integrity and quality of an image. However, the results have much wider 
implications in the generic behaviour of fluid droplets on heterogeneous sur- 
faces. This is of particular interest since patterned surfaces are used in wide 



1 



diversity of areas such as electronic devices, biologically active substrates or 
semiconductor nanostructures. 

Previous work on the behaviour of fluids on chemically patterned sub- 
strates has predominantly concentrated on the equilibrium shape of the 
dropsP In particular there has been extensive theoretical and experimental 
work on the extent to which the line tension affects droplet shape] EElSIZI g ev _ 
eral authors have considered the equilibrium configuration of a fluid on one or 
two completely wetting stripes in a non-wetting 

there is recent work describing how such a fluid moves along such a stripe.^ 
By contrast, this paper describes the behaviour of fluid droplets on chemi- 
cally patterned surfaces when the fluid is jetted onto the surface (non-zero 
impact velocity). This offers new perspectives regarding the achievement of 
structures impossible to obtain by simple droplet deposition. 

The experimental results are compared to numerical simulations achieved 
by a lattice Boltzmann solution of the equations of motion for a one-component, 
two-phase fluid. Lattice Boltzmann models are a class of numerical tech- 
niques ideally suited to probing the behaviour of fluids on mesoscopic length 
scales.^ Several lattice Boltzmann algorithms for a liquid-gas system have 
been reported in the literature! 18 * 19 * 20 * They solve the Navier-Stokes equa- 
tions of fluid flow but also input thermodynamic information, typically either 
as a free energy or as effective microscopic interactions. They have proved 
successful in modelling such diverse problems as fluid flows in complex geome- 
tries,^ two-phase modelsp^^ hydrodynamic phase ordering^ and sediment 
transport in a fluidP^ 

We consider a one- component, two-phase fluid and use the free energy 
model originally described by Swift et al. l8: with a correction to ensure 
Galilean invarianceP^ The advantage of this approach for the wetting prob- 
lem is that it allows us to tune equilibrium thermodynamic properties such as 
the surface tension or static contact angle to agree with analytic predictions. 
Thus it is rather easy to control the wetting properties of the substrate. 

Choosing fluid parameters such as viscosity, surface tension and contact 
angles to match the experimental values we obtain a close quantitative match 
to the different droplet shapes thus helping to understand the mechanisms for 
their formation. The results show the importance of dynamic and metasta- 
bility effects in determining the shapes formed by small droplets on het- 
erogeneous substrates. In particular we demonstrate that the final droplet 
configuration depends on the droplet size relative to that of the stripes, the 
initial point of impact and the incident droplet velocity. 
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2 Experimental Section 



2.1 Fluid 

The solventless UV cure black ink-jet ink used (acrylate monomer) has a 
viscosity of 25.0 mPa.s and a surface tension of 24.4 mN.m" 1 . The droplets 
were jetted using a 256 nozzle industrial inkjet printhead held 1 mm from the 
substrate. The print pattern was configured so that the droplets impacted 
on the surface with an average separation between each drop of 280 /im, at a 
velocity of 8 m.s -1 . The ink droplets were cured 1.6 s after ejection from the 
print head using a standard mercury H UVA lamp system. The radius of the 
droplets before impact was chosen to be 22 /zm, and the typical spreading 
time of the droplets is of order milliseconds. 

2.2 Preparation of the substrates 

Surfaces are produced with areas of different wettabilities using standard mi- 
crocontact printing techniques. This is nowadays a well established method, 
and abundant literature can be found on the subject.^ The chemical pat- 
terns are created on gold coated silicon (001) wafers to make certain that the 
substrates are molecularly flat, thereby ensuring any effects observed in the 
droplet behaviour resulted from chemical rather topographic (or a combina- 
tion of both) effects. Prior to gold evaporation, the wafers are coated with 
50 nm chromium to ensure stability of the gold (100 nm). The samples are 
washed with ethanol and dried with a stream of N 2 before patterning. 

The stripes are created using solutions of 4 mmol of either methyl- or 
carboxyl- terminated octadecyl thiols, which produced lyophobic (-CH3) and 
lyophilic (-COOH) monolayer regions on the Au surface. To achieve this, 
a polydimethylsiloxane stamp is moulded on a master having the desired 
pattern, and allowed to cure for one week at room temperature. Then a 
known quantity of hexadecanethiol solution (4 mmol in hexane) is poured 
onto the elastomeric stamp. After 10 seconds and drying with a stream of 
N 2 , the stamp is applied onto the surface of the wafer. 

This produces a lyophobic self assembled monolayer (SAM) with the de- 
sired pattern. The sample is then immerged one hour in a mercapto unde- 
canoic acid / ethanol solution (4mmol) in order to complete the patterning. 

The corresponding contact angles for the commercially available black 
UV cure jet-ink used in these experiments are respectively 64° and 5° on the 
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Table 1: Widths of the chemically patterned stripes used in this study, as 
determined by scanning electron microscopy. The sample numbering corre- 
sponds to that used in Figure [T] 

Sample lyophilic (/iin) lyophobic (/xm) 



SI 


5 


56 


S2 


13 


65 


S3 


19 


35 


S4 


26 


47 


S5 


23 


38 


S6 


50 


50 


S7 


47 


31 


S8 


83 


75 



lyophobic and lyophilic areas. The contact angles are determined by optical 
measurement of the equilibrium shape of equivalent sized droplets on surfaces 
homogeneously produced by solution casting the relevant thiol. The widths 
of the stripes created by the microcontact printing process, and measured 
using scanning electron microscopy, are listed in Table ^ 

3 The lattice Boltzmann model 

The lattice Boltzmann approach solves the Navier-Stokes equations by fol- 
lowing the evolution of partial distribution functions fi on a regular, d- 
dimensional lattice formed of sites r. The label % denotes velocity directions 
and runs between and z. DdQz + 1 is a standard lattice topology classi- 
fication. The -D3Q15 lattice topology we use here has the following velocity 
vectors v*: (0,0,0), (±1,±1,±1), (±1,0,0), (0,±1,0), (0,0, ±1) in lattice 
units. 

The lattice Boltzmann dynamics are given by 

/ 4 (r + Atv„ t + At) = Mr, t) + - (/f(r, t) - /,(r, t)) (1) 

r 

where At is the time step of the simulation, r the relaxation time and f^ q 
the equilibrium distribution function which is a function of the density n = 
Y%=o ft an d the fluid velocity u defined through the relation nu = Y%=o fi v i- 
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The relaxation time tunes the kinematic viscosity as^ 

where Ar is the lattice spacing and C 2 and C4 are coefficients related to the 
topology of the lattice. These are equal to 3 and 1 respectively when one 
considers a D3Q15 lattice (see^for more details). 

It can be showrP^ that equation (JTJ) reproduces the Navier-Stokes equa- 
tions of a non-ideal gas if the local equilibrium functions are chosen as 

fl q = A a + B a u a v ia + C a u 2 + D a u a U{3V ia Vi(3 + G aa pV ia Vip, i > 0, 

fo q = n-E/r (3) 

i=i 

where Einstein notation is understood for the Cartesian labels a and (3 (i.e. 
ViaU a = J2 a v iaU a ) and where a labels velocities of different magnitude. A 
possible choice of the coefficients is^ 

A a = (p b - ^(d a n) 2 - nnd aa n + vu a d a nj , 



C 2 



2c 2 ' 2c 4 



^177 = (n^n) 2 + 2vu 1 d 1 n s j , G 277 = 0, 

G 2l s = ( K (9 7 n)(d s n) + u(u 7 d s n + u s d 7 n)) (4) 

where W\ = 1/3, w 2 = 1/24, c = At/ At, k is a parameter related to the 
surface tension and pb = p c (y p + l) 2 (3z/ 2 — 2u p + 1 — 2f3r p ) is the pressure in 
the bulk where v p = (n — n c )/n c , t p = (T c — T)/T c and p c = 1/8, n c = 3.5 
and T c = 4/7 are the critical pressure, density and temperature respectively 
and (5 is a constant typically equal to 0.1. 

The derivatives in the direction normal to the substrate are handled in 
such a way that the wetting properties of the substrate can be controlled. A 
boundary condition can be established using the Cahn modelP^ He proposed 
adding an additional surface free energy ^ c (n s ) = 4>o — 4>\n s + ■ • • at the solid 
surface where n s is the density at the surface. Neglecting the second order 
terms in \I/ C (n) and minimizing ^ + \l/ c (where is the free energy in the 
bulk), a boundary condition valid at z = emerges 

d z n = - — . (5) 

K 
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Equation (jSJ is imposed on the substrate sites to implement the Cahn 
model in the lattice Boltzmann approach. Details are given inP^ 

The Cahn model can be used to relate 0i to 6 the contact angle defined 
as the angle between the tangent plane to the droplet and the substrate^ 

0i = 2(3r p ^2p c K sign(# - |)^cos| - cos (6) 

where a = cos _1 (sin 2 6) and the function sign returns the sign of its argu- 
ment. 

We impose a no-slip boundary condition on the velocity. Because the full 
dynamics takes place on the boundary the usual bounce-back condition must 
be extended to ensure mass conservation (see^for a wider discussion). This 
is done by a suitable choice of the rest field, /o, to correctly balance the mass 
of the system. 

This model reproduces Young's law and the expected dependence of the 
droplet behaviour on viscosity and surface tension.^ 

The following lattice Boltzmann parameters are set. The initial droplet 
radius Ro = 30 lattice sites. The droplet is initialised with a vertical velocity 
equal to U = 0.02. The lattice geometry is L x x L y x L z where L x and L y 
are chosen large enough to not affect the behaviour of the droplet and L z = 
40. The relaxation time r = 0.63. The surface tension related parameter 
k = 0.0012. The temperature T = 0.4 which leads to two phases of density 
ni = 4.128 and n g = 2.913. The simulations are run for 400 000 iterations. 

Simulation and physical parameters are related as usual by choosing a 
length scale L Q , a time scale T and a mass scale M . A simulation parame- 
ter with dimensions [L} ni [T] n2 [M] na is multiplied by L^T^Mq 3 to give the 
physical value. L = 1.5 • 10~ 6 m, T = 3.7 • 10~ 9 s and M = 7.9 • 10~ 16 ^ 
were chosen to give physically realistic droplet diameter, viscosity and surface 
tension. 

4 Results and discussion 

A top view of the final droplet shapes on the different patterned surfaces is 
shown in Figure ^ To our knowledge, it is the first time that the evolution 
of droplet shape has been studied over such a large range of heterogeneities. 

When the width of the stripes is much smaller than the droplet radius 
(5*1) the contact line is almost circular with only small deviations near the 
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Figure 1: Scanning electron micrographs of inkjet droplets on patterned 
surfaces. Lyophilic and lyophobic stripes appear dark and pale respectively."*! 
The stripe widths for the different saiLples are listed in Tabled 



stripe boundaries as previously reported in the literature)^ 

For wider stripes, where the dimensions are now approaching the droplet 

diameter, the fluid remains confined to a lyophilic region (S6, S8) and the 

final shape is therefore highly elongated. 

It should also be noted that in samples S6 and S8, the regular droplet 

pattern array is disrupted as droplets can be displaced to the neighbouring 

lyophilic region. When this occurs, a small portion of the initial droplet 

volume can remain on a neighbouring lyophilic stripe. 

In the intermediate cases, for droplets impinging on the stripes S3, S4, S5 

and 57, two different characteristic droplet shapes are observed: "lozenge" 

or "butterfly". 

The reason for the formation of these shapes can be understood by con- 
sidering a lattice Boltzmann modelling of the droplet dynamics. 

Figure Efa) shows the droplet contact lines for intermediate and equilib- 
rium stages of the droplet wetting as a function of the initial impact position 
of the jetted fluid. The modelled substrate is defined so that the surface 
heterogeneities are equivalent with the sample 5*4. The droplet at the point 
of impact is assumed to be spherical with an imposed velocity equal to that 
used experimentally (8 m.s -1 ). It is quite clear that the shapes obtained by 
such calculations are entirely consistent with those obtained experimentally. 
The very strong similarities are more clearly demonstrated in Figure Efb) 
where the experimental and numerical simulations are overlaid for the two 
equilibrium droplet shapes. The good quantitative agreement between the 
simulations and the experimental results is extremely pleasing given that all 
the parameters in the simulation are fixed by the physical conditions dictated 
by the experiments. Small differences are not surprising because of the un- 
certainty in the experimental values of the transport coefficients and surface 
energies, and the possibility of contact line pinning on real substrates. 

Note in particular that, just as in the experiments, two different final 
droplet shapes are obtained in the numerical modelling for the substrate 
geometry S4. 

The simulations allow us to follow the dynamical pathway by which the 
final equilibrium states are reached. Consider case 1 (see Figure Ofa)) where 
the initial contact point of the droplet is in the centre of a lyophilic stripe. 
The drop initially evolves symmetrically parallel and normal to the stripe. 
Once the boundary of the lyophobic stripe is reached horizontally, the fluid 
moves more quickly along the lyophilic than along the lyophobic part of 
the substrate. Note that there is an overshoot of the fluid in the direction 
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perpendicular to the stripes (the maximum extent of the contact line exceeds 
that of the equilibrium droplet). This occurs as the droplet initially spreads 
due to its impact velocity, and relaxes to its equilibrium shape determined 
by the wetting behaviour. The net effect is an equilibrium shape which has 
a distinctive lozenge shape centred over a lyophilic stripe. 

Now consider case 2 where the droplet strikes the substrate in the middle 
of a lyophobic stripe. Because the initial diameter is approximately the same 
as that of the lyophobic stripes, the droplet spreads directly onto adjacent 
lyophilic parts of the substrate giving the characteristic butterfly shape which 
is symmetrically located over two lyophilic and one lyphobic stripes. 

When the droplet impacts on a location other than the middle of a stripe, 
the symmetry of the evolution parallel to the stripes is lost. Despite this 
initial asymmetry the final state is always symmetric and produces either 
the characteristic lozenge (for the initial positions 3, 5 and 6) or butterfly 
(for the initial position 4) patterns. These states correspond to long lived 
metastable or stable equilibria.^ The numerical results give a free energy 
for the butterfly-shaped droplets which is ~ 1% of the lozenge free energy. 

Figure |2fa) also allows us to predict that jetting droplets onto substrate 
54 should lead to approximately equal numbers of lozenge and butterfly 
shapes. Numerical results show that a droplet spreads into a lozenge if it hits 
the surface on a lyophilic stripe or its close surroundings (up to about one 
eight of the width of the neighbouring lyophobic stripes). This corresponds 
to about half the substrate area and therefore one expects about one half of 
the droplets to take each of the final shapes consistent with the experimental 
results for 54, as shown in 

Figure H Thus far we have shown that the parameters affecting the final 
shape of a drop impacting on a striped substrate are the relative sizes of the 
drop and the stripes and the initial point of impact. 

Clearly, the surface energies of the different stripes will have an important 
effect on the droplet spreading behaviour, U: and given the results described 
above the initial drop velocity may also be important. 

To test the latter, numerical simulations comparing the time evolution of a 
drop with and without impact velocity have been undertaken (see Figure EJ). 
The final state is indeed different, as the drop with no impact velocity is 
unable to reach, and hence take advantage of the wetting possibilities of the 
neighbouring lyophilic regions of the substrate. For a given average distance 
between heterogeneities, the droplet velocity appears to be an important 
parameter regarding to the number of surface heterogeneities encountered 
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Figure 2: The modelled substrate is defined so that the surface hetero- 
geneities are equivalent with the sample S4. (a) Numerical simulation of 
droplets hitting the surface at various impact points indicated by encircled 
crosses. For each droplet the bold and faint lines represent the extent of the 
droplet at equilibrium and at intermediate times, respectively. The lyophobic 
and lyophilic regions are shaded to be consistent with Figure [TJ (b) Direct 
comparison between experimental (S4, dashed lines) and numerical (solid 
lines) equilibrium contact lines. The length scales in both plots are reported 
in micrometers. 




10 



(a) (b) 




Figure 3: Time evolution (top to bottom) of the droplet shape from numerical 
simulations of the equations of motion, (a) Droplet spreading with impact 
velocity 8 ms _1 . (b) No impact velocity. The initial point of impact is the 
same for both simulations, tf is the final simulation time. 
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during spreading. 



5 Conclusion 

We have presented experimental and numerical results investigating the be- 
haviour of micron-scale droplets on chemically patterned substrates. When 
the drop radius is of the same order as the stripe width the final droplet shape 
is determined by the dynamic evolution of the drop and is very sensitive to 
the initial droplet position and velocity. The final state may be metastable 
and is not determined just by a minimisation of the free energy. 

We have shown that it is possible to provide a close quantitative corre- 
spondence between numerical solutions of the hydrodynamic equations of 
motion describing the spreading and the experimental results. This has 
proved invaluable in fully understanding the data and in predicting droplet 
behaviour for parameter values not available experimentally. 

These results underline the difficulties inherent in controlling the details 
of patterns formed using ink-jet printing (where the underlying substrate is 
likely to have both chemical and topological heterogeneities) and the subtle 
effects of the surface wetting properties on the behaviour of liquids on pat- 
terned substrates. The present work will therefore be extended to topological 
surface features and complex fluids. 
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